Five-year Monitorings of TeV Blazars with ASCA and RXTE 
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Abstract 

We study the temporal/spectral variability of two cxtragalactic TeV sources, Mrk 421 and Mrk 501, 
based on 5-year observations with the ASCA and RXTE satellites. We found that the peak of the syn- 
chrotron emission exists just in the X-ray band and its position shifted from lower to higher energy when 
the source became brighter. Relationship between the peak energy and peak luminosity showed quite 
different behavior in the two sources; Mrk 421 showed very little change in the peak position (0.5—2 keV), 
while Mrk 501 revealed the largest shift ever observed in blazars (1 — 100 keV). We analyze these X-ray 
data with the flux changes in TeV band, which are obtained from 35 truly simultateneous observations. 
Very different spectral evolution of both objects indicates some differences in the electron acceleration 
mechanism at work during the flares. We argue that the flux variability of Mrk 421 is associated with 
an increase in the number of electrons, while the flare of Mrk 501 is mostly due to the large changes in 
maximum energy of electrons. We also discuss the characteristic temporal variability of TeV sources, and 
implications for the X-ray emitting site in the relativistic jet. 

Key WORDS: blazars: time variability — spectrum: evolution 



1. Introduction 

Blazars, a subclass of AGNs, have outstanding proper- 
ties in several aspects. They exhibit the most rapid and 
the largest amplitude variations of all AGNs. Recent 
observations with the EGRET instrument on-board the 
Compton Gamma-Ray Observatory {CGRO) reveal that 
more than 60 AGNs are bright 7-ray emitters (e.g.. Hart- 
man ct al. 1999). Remarkably, most of AGNs detected 
by EGRET were classified as blazars. Observations with 
ground-based Cherenkov telescopes further confirmed 7- 
ray emission extending up to TeV energies for a number 
of nearby blazars. 

Overall spectra of blazars (plotted as v F„) have at 
least two pronounced continuum components: one peak- 
ing between IR and X-rays, and another in the 7-ray 
regime (e.g., von Montigny et al. 1995). The strong po- 
larization observed in the radio and optical bands (Angel 
& Stockman 1980) implies that the lower energy compo- 
nent is produced by synchrotron radiation of relativis- 
tic electrons in magnetic fields, while inverse-Compton 
scattering by the same electrons is dominant process 



responsible for the high energy 7-ray emission (Ulrich, 
Maraschi, & Urry 1997). This non-thermal radiation is 
thought to be emitted in a relativistic jet pointing close 
to our line of sight (e.g., Urry & Padovani 1995). 

Multi-frequency spectra of blazars indicate their vari- 
ety. There are indications that particles are accelerated 
more efficiently in blazars having lower luminosities (e.g., 
Ghisellini et al. 1998). In fact, the luminosity of TeV 
blazars is the lowest end of EGRET detected blazars, 
while the emission peaks are located at the highest fre- 
quencies. In low luminosity objects, energy loss by radi- 
ation would not be effective because there are less ambi- 
ent photons to be Compton scattered to higher energies. 
Thus the TeV blazars gives an important opportunity to 
study the particle acceleration in blazar jets. 

Mrk 421 and Mrk 501 are the proto-type of TeV 7-ray 
emitting blazars. The first multi-frequency campaign of 
Mrk 421 including TeV energies was conducted in 1994 
(Macomb et al. 1995; Takahashi et al. 1996). The con- 
temporaneous observations implied correlated variability 
between the keV X-ray and TeV 7-ray emission, while 
the GeV fiux and the radio to UV fluxes showed less vari- 



ability. The results from this campaign are important 

because it suggested the possibility that a single elec- 
tron population is responsible for both the X-rays and 
TeV 7-rays, qualitatively agreeing with the Synchrotron 
self-Compton (SSC) scenario. However, the fact that the 
X-ray and TeV observations were separated by one day, 
could introduce uncertainties to draw concrete pictures. 
Simultaneous, uninterrupted observations are crucial to 
understand the blazar phenomenon. 

In this paper we intensively study Mrk 421 and Mrk 
501, with an emphasis on the observations with the X-ray 
satellites ASCA and RXTE . Unlike most of the previ- 
ous work which focused on the photon spectral properties 
and/or classification of blazars, we also study the rapid 
time variability and spectral evolution in blazars. The 
time variability gives us important constraints on the 
dynamics operating in blazar jets, such as the accelera- 
tion, radiative cooling and particle escape. Taking the 
VLBI observations of super luminal motion into account, 
we discuss the internal jet structure and emission site in 
the jet. 

2. X-ray Observations and Data Analysis 
2.1. Observations 

Mrk 421 and Mrk 501 were observed a number of times 
with the X-ray satellites ASCA and/or RXTE . ASCA 
observed Mrk 421 five times with a net exposure of 
546 ksec between 1993 and 1998. In the 1998 obser- 
vation, the source was in a very active state and was 
detected at its highest-ever level (Takahashi et al. 2000). 
RXTE ohseived Mrk 501 more than 100 times with a net 
exposure of 700 ksec between 1996 and 1998. Mrk 501 
was in a historical high-state in 1997 (Catanesc ct al. 
1997; Plan et al. 1998; Lamer & Wagner 1999). Multi- 
wavelength campaigns, including a number of ToO (Tar- 
get of Oppotunity) observations, were conducted during 
this high state. Observation logs are given in Table 1. 



Table 1. Observation log of Mrk 421 and Mrk 501 









Exp. 


ID 


Source 


Satellite 


Obs. Time(UT) 


(ksec) 




Mrk 421 


ASCA 


93.05.10 - 05.11 


43 


(a) 




ASCA 


94.05.16 - 05.17 


39 


(b) 




ASCA 


95.04.25 - 05.08 


91 


(c) 




ASCA 


97.04.29 - 05.06 


70 


(d) 




ASCA 


98.04.23 - 04.30 


280 


(e) 


Mrk 501 


RXTE 


97.04.03 - 04.16 


36 


(f) 




RXTE 


97.05.02 - 05.15 


51 


(g) 




RXTE 


97.07.11 - 07.16 


38 


(h) 




RXTE 


98.02.25 - 10.02 


568 


(i) 



2.2. X-ray Time Variability 

Figure 1 shows the long-term variation of fluxes, with 
both ASCA and RXTE observations plotted. ASCA ob- 
servations of Mrk 421 spanned more than 5 years (from 
1993 to 1998), and show that the source exhibits vari- 
ability by more than an order of magnitude. Blow-ups 
of the light curves taken in 1995 and 1997 are given in 
the lower panel (Figure 1 (c),(d)). For Mrk 501, we plot 
the RXTE data because the observations were conducted 
much more frequently than the ASCA observations. The 
RXTE observations spanned more than 3 years and sig- 
nificant flux changes are clearly detected. Blow-ups of 
light curves are given in Figure 1 (f) — (h). 

To cvahiatc variability features in blazars, we exam- 
ine the use of a numerical technique called the structure 
function (hereafter, SF). While in theory the SF is com- 
pletely equivalent to traditional Fourier analysis meth- 
ods, it has several significant advantages. Firstly, it is 
much easier to calculate. Secondly, the SF is less af- 
fected by gaps in the light curves. The definitions of SFs 
and their properties are given by Simonetti et al. (1985). 
The first order SF is defined as 

SF{r) = j^Y.^ait)-ait + r)r, (1) 

where a{t) is a point of the time series (light curves) {a} 
and the summation is made over all pairs separated in 
time by t. N is the number of such pairs. 

The structure functions calculated from the light 
curves (Figure 1) are given in Figure 2. The SFs are nor- 
malized by the square of the mean fluxes, and are binned 
at logarithmically equal intervals. Both the SFs show a 
rapid increase up to r/day ~ 1, then gradually flatten to 
the observed longest time-scale of r/day > 1000. Fluctu- 
ations at large r (r/day > 10) are due to the extremely 
sparse sampling of data. In order to demonstrate the un- 
certainties caused by such sparse sampling, and to firmly 
establish the reality of the "roll-over" , we simulated the 
long-term light curves following the forward method de- 
scribed in lyomoto (1999). We found that (1) the PSD 
of the TeV sources have at least one roll-over at 10~^ Hz 
< /br < 10-^ Hz (1 < r/day < 10), and (2) the PSD 
changes its slope from oc J-i~-2 ^ j^^-j ^ j-2 3 
(/ > ./br) around the roll-over (Kataoka et al. 2001a). 

2.3. X-ray Spectral Evolution 

We next performed model fitting to evaluate the X-ray 
photon spectra of Mrk 421 and Mrk 501. In general, 
the photon spectra of TeV sources show convex shape, 
where the model with a single power law functions plus 
absorption arising from neutral material could not fit the 
observed spectra. We thus use another model which is 
expected to fit the spectrum 

^ = NoyieM-N^ME))E-^^^M-EclEU2) 
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Fig. 1. Long-term flux variation of Mrk 421 (left) and Mrk 501 (right) 



where T is the photon index and Ec is the cutoff energy, 
respectively. Using this relatively simple model, one can 
differentiate the function analytically. The peak energy 
Ep in v Fi, are given E^ = We cannot determine 

Ep for r < 2, because in such cases, the spectra are 
monotonously rising in the total energy band. 

We divided the observations into 5 kscc cxporsure seg- 
ments and fitted the photon spectra for individual seg- 
ments. In this case, statistically acceptable fits were ob- 
tained. The distribution of peak luminosity Lp versus 
peak energy Ep is summarized in Figure 3. The peak 
luminosity is simply calculated from the luminosity at 
E = Ep. Mrk 421 shows very modest shifts in the peak 
position. A clear correlation between the peak energy 
Ep and the peak luminosity Lp was detected for the first 
time for both Mrk 421 and Mrk 501: Ep oc L^-^ for Mrk 
421 and Ep oc L^-^ for the case of Mrk 501. Extrapo- 
lating the relation between Ep and Lp, we expect that 
synchrotron peak reaches to ~ 100 keV in the high state 
of Mrk 501, which is the largest shift ever observed in 
blazars (Plan et al. 1998; Kataoka et al. 1999). 

3. Multi-frequency Analysis 

In previous sections, wc have investigated the rapid vari- 
ability and spectral evolution of TeV blazars, highlight- 
ing the X-ray observations by ASCA and RXTE . How- 
ever, the overall spectral energy distribution of blazars 
generally ranges over a very wide range - from radio 
to TeV 7-ray bands - and such spectra are one of the 
most important features allowing us to understand these 
sources. To reveal the multi-frequency properties of TeV 
blazars, we have to carry out multi-frequency monitoring 
campaigns of the source in various states of activity. 



Results from previous campaigns strongly suggest cor- 
related variations in X-ray and TeV 7-ray flirxes, while 
variability is much less pronounced in other energy 
bands. However, a problem remains that most of the 
data are taken either non-simultaneously or are very 
sparsely sampled (e.g., Macomb et al. 1995; Buckley 
et al. 1996). The discussion based on those quasi- 
simultaneous data may be incomplete. We need exactly 
simultaneous monitoring, especially in the X-ray and 
TeV energy bands, to correctly understand the sources. 
In the following, we summarize the results from truly 
simultaneous campaigns of TeV blazars conducted from 
1996 to 1998. 

Evolution of multi- frequency spectrum of Mrk 421 is 
shown in Figure 4:{left). Open circles come from Ma- 
comb et al (1995), while other symbols represent our 
new results taken exactly simultaneously in X-ray (also 
EUVE) and TeV 7-ray bands. Several important fea- 
tures are seen in the figure. First, a clear correlation 
can be seen between the keV X-ray flux and the TeV 
7-ray fluxes. The amplitude of TeV flux variation is 
less or almost comparable with that in the X-ray energy 
bands. Detailed analysis showed that the correlation is 
expressed as [TeV 7-ray flux] oc [X-ray flux]'^-^^='=°--^^. Sec- 
ond, the slope of X-ray photon spectra (Fx-ray ~ 3) are 
steep, and very similar to those in the TeV energy band 
(FxeV ~ 3; Aharonian et al. 1999b). Third, although 
X-ray flux changed dramatically in various seasons, only 
small changes are implied in the synchrotron peak posi- 
tion, as wc have quantitatively discussed in § 2.3. 

Evolution of the multi- frequency spectrum of Mrk 501 
is shown in Figure 4 (right). Open circles are data from 
public archive (taken from Kataoka et al. 1999), while 
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Fig. 3. Distribution of peak luminosity Lp versus peak energy Ep 



other symbols represent our new results derived from 
simultaneous monitoring in the X-ray and TeV 7-ray 
bands. For the data obtained in March 1996, we also 
plot the EGRET and optical data. Compared to the 
results for Mrk 421, some significant differences are im- 
plied in the multiband specta. First, the changes in X- 
ray flux are accompanied by large shift in the position of 
the synchrotron peak. Second, the slope of X-ray spec- 



tra (F) varies widely, ranging from Fx- 



1.7 to 2.5, 



while the photon index in the TeV energy band is almost 
unchanged; FxeV ~ 2.5 (e.g, Aharonian et al. 1999a). 
Third, amplitude of TeV flux variation is much larger 
than that in the X-ray energy band. The correlation be- 
tween two energy bands is expressed as [TeV 7-ray flux] 
oc [X-ray flux]i-96±o.07_ 



4. Discussion: what can we learn? 

As first pointed out by Kataoka (2000), the little power 
of rapid variability (< 1 day) in TeV sources provides 
important clues to the X-ray emitting site in the jet. 
The characteristic time-scale of each flare event ivar > 1 
day should reflect the size of the emission region, which 
we infer to be > 10^^ (F/IO) cm in the source co-moving 
frame, if emitting blobs are approaching with Lorentz 
factors F (F ~ 10; Vermeulen & Cohen 1994). This 
range of Lorentz factors is independently inferred from 
constraints that can be derived from the spectral shape 
and variability of TeV blazars (e.g., Tavecchio, Maraschi 
and Ghisellini 1998; Kataoka et al. 1999). 

If the jet is coUimated to within a cone of constant 
opening angle 6 ~ 1/F and the line-of-sight extent of 
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Fig. 4. Evolutions of multi-frequency spectra : Mrk 421 (left) and Mrk 501 {right). Dotted line represents the best fit SSC model for the 
quiescent/flare state. Figure from Kataoka et al. (2001b). 



shock is comparable with the angular extent of the jet, 
one expects that the X-ray emission site is located at 
distances D > 10^^ (F/IO)^ cm from the base of the jet. 
Only little variability shorter than fvar strongly suggests 
that no significant X-ray emission can occur in regions 
closer than this to the black hole. The relativistic elec- 
trons responsible for the X-ray emission are most likely 
accelerated and injected at shock fronts occurring in the 
jet (e.g., Inoue & Takahara 1996). The lack of short term 
variability may then imply that shocks are nearly absent 
until distances of I? > lO^"^ (r/lO)^ cm. Further discus- 
sion of the characteristic timescale is given in Kataoka et 
al. (2001a), which also suggests a possible link between 
the jet and the central engine. 

In § 2.3, we found that the position of the synchrotron 
peak shifts from lower to higher energy when the source 
becomes brighter. The difference of spectral evolution in 
Mrk 421 and Mrk 501 implies that quite different mecha- 
nisms were at work when the sources went into the flaring 
states. Since the peak luminosity is proportional to the 
number of photons at peak (nph(£^p)) multiplied by the 
peak energy, a simple relation can be found: 

Lp (X _Bpnph(£^p) cx 7pno(7p), (3) 

where 7p is the Lorentz factor of electron which emits 
photon of energy Ep and ne(7p) is the number of elec- 
trons at 7p. We thus find the relations, 

ne(7p) cx 7^-' (Mrk421), (4) 
ne(7p) (X 7p (MrkSOl). (5) 

This implies that for the case of Mrk 421, increase of 7p 
by a factor of 2 requires more than factor 10 increase in 
number of electrons at peak energy. On the other hand. 



the same amount of increase of 7p requires decreases 
in number of electrons by a factor of 2 for the case of 
Mrk 501. Such variations in spectral behaviors may be 
associated with the difference of physical conditions in 
relativistic jets, which is discussed below. Full descrip- 
tions are given in Kataoka (2000) and Kataoka et al. 
(2001b). 

Assume that a substantial amount of gas (e.g., in a 

form of clouds) is distributed in the jet. In these clouds, 
density of low-energy electrons is enhanced as compared 
to the ambient, but other physical quantities such as 
magnetic field strength arc unchanged. Those clouds 
essentially provide a plentiful source of low-energy elec- 
trons for the shock front. When a shock front overruns 
one of such clouds, fresh electrons are successively in- 
jected and assumed to undergo continuous acceleration 
by repeatedly crossing and recrossing the shock front, 
as well as simultaneously cooling by synchrotron radi- 
ation (e.g.. Kirk, Rieger & Mastichiadis 1998). In this 
scenario, number of electrons increases significantly, but 
only small changes are implied for the maximum Lorentz 
factor as were observed in Mrk 421. 

On the other hand, when clouds are absent or very 
sparsely distributed in the jet, fiares may be produced 
in several diffent manners. For example, if the shock 
overruns the enhanced tangled-magnetic field region, this 
may cause changes in acceleration time of electrons, 
and hence increase the maximum Lorentz factor. Im- 
portantly, regardless of detailed models for flaring be- 
haviour, total number of electrons is conserved in this 
case. During the flare, acceleration can be assumed to 
be more efficient than radiative cooling, thus the present 
electron population as a whole will be accelerated to 
higher energies, but no additional electrons are supplied 



into the shock. 

In latter case, numbers of electrons at the peak (tIc 
(7p)) ^ill decrease, reflecting the power-law shape of an 
electron population. If the differential number density 
of electrons is expressed as iVo oc -f~^ (standard shock), 
electron number decreases as nd^p) oc 7p-/Ve(7p) oc 7~^. 
Importantly, this relation is very close to the case we 
have observed in Mrk 501. 

Present discussion based on the X-ray spectral evo- 
lution suggests very important implications for internal 
jet structures. Only a small shift of synchrotron peak 
observed in Mrk 421 may be associated with electron 
clouds filling the jet, while the jet of Mrk 501 seems 
to be relatively empty. During the flare of Mrk 421, ki- 
netic power of the shock is equally distributed to large 
number of low-energy electrons newly injected into the 
shock, thus increasing the number of high energy elec- 
trons. Large shifts of synchrotron peak observed in Mrk 
501, on the other hand, is possible only when the internal 
jet is rather sparse and transparent to the shock propa- 
gation. Kinetic power of the shock is spent to increase 
the energies of individual electrons and hence number- 
conservative. 

It may thus be worthwhile to compare our X-ray im- 
plications to the VLBI results. Accurate measurements 
of changes in the parsec-scale jet structure imaged with 
VLBI provide constraints on the jet kinematics and ge- 
ometry. When combined with estimates of the Doppler 
beaming factor (determined, for example, from the X-ray 
time variability), the apparent motion of the jet compo- 
nents can be used to constrain the Lorentz factor of the 
jet and the angle of the jet to the linc-of-sight. Although 
VLBI observations do not yet have an enough resolution 
to image the region production of the X-ray/7-rays (~ 
0.01 pc), they provide the highest resolution structual 
information available, and can image the region immedi- 
ately downstream. 

Mrk 421 and Mrk 501 are observed in a space VLBI 
pro j ect using the HALCA sa.tdlite and 12 ground sta- 



421 and Kataoka et al. (1999) derived (5 > 6 for Mrk 501 
using the X-ray/TeV 7-ray time variability. Marscher 
(1999) disscussed that rather low superluminal apparent 
speeds and lackluster variability properties of the radio 
jets evidences that the bulk flow of the jets decelerates 
from X-ray/TeV 7-ray emitting section (~ 0.01 pc) to 
the radio emitting region (~ 1 pc). The weak subparsec- 
and parsec-scale jets of those objects are readily under- 
stood as the consequence of heavy energy and monmen- 
tum loss in the upstream of the jet where most of the 
energy and momtentum of the relativistic electrons are 
transfered to the radiation in the X-ray and TeV 7-rays. 

Remarkably, very different manners of X-ray spectral 
evolutions of Mrk 421 and Mrk 501 presented in this pa- 
per is exactly consistent with those VLBI observations. 
Our X-ray observations predict that relativisitc outflows 
of Mrk 421 will be decelerated faster than that for Mrk 
501, because the jet of Mrk 421 is filled with low-energy 
materials and kinetic energy of outflows are more effi- 
ciently dissipated during the propagation. The absence 
of superluminal motion of Mrk 421 in radio bands thus 
indicates that the high energy outfiow has been suffi- 
ciently decelerated when it reaches to more distant, radio 
emitting region. 
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